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ABSTRACT 

1 1 

The solution of the one-dimensional time-independent Schrodinger 
equation in which the energy minus the potential varies as the n-th 
power of the distance is obtained from proper linear combinations of 
Bessel functions,, of order (n + 2) \ The linear combinations whieh^we 
call^'generalized Airy functions" (2jL'^(x ) ■- -• and i3i^(x) , reduce to the 
usual Airy functions Ai(x) and Bi(x) when n - 1 and have the same 
type of simple asymptotic behavior. Expressions for the generalized 
Airy functions which can be evaluated by the method of generalized 
Gaussian quadrature are obtained. 
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GENERALIZED AIRY FUNCTIONS FOR USE IN ONE-DIMENSIONAL QUANTUM MECHANICAL 

PROBLEMS 


For one-dimensional quantum mechanical systems, the usual WKB 
approximation gives simple expressions for the wave function in both 
the classical and the non-classical regions of configuration space. 
Connection formulas relate the two expressions on either side of a 
classical turning point where the WKB wavefunctions become infinite. 

The standard textbook^ example of the connection formulas in the semi- 
classical approach assumes that the potential is a linear function of 
x in the neighborhood of the turning point. The exact solution in 
this case is given by the Airy functions Ai(x) and Bi(x) which are 
simply related to Bessel functions of order 1/3. An analogous treatment 
for the Schrodinger equation 

2 

d w(x) , 2 n , x ,, N 

x — + c x w (x) = 0 , (1) 

dx 


where c^x 11 = Q^(x) = ^ [E - V(x)] has a zero of order n is easily 

-fi 

obtained but does not seem to be readily available. In this note explicit 

2 

expressions for the solution of Eq. (1) for which Q has a zero of 

arbitrary order are obtained. We call these solutions ’’generalized 

Airy functions" (Xi^Cx) and H3i (x) since they are standardized to 

2 

give the Airy functions for n = 1 . Gordon has given a method for 
calculating the Airy functions by means of generalized Gaussian quadra- 
ture. This method is extended here for the calculation of the generalized 
Airy functions . 
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In the following development it is shown that the generalized Airy 
functions are linear combinations of solutions w^ + (x) to Eq. (1), 


(Xi v (x) = [w v + (x) + w v (x)] 

l6i v (x) = k 2 [w v (x) - w v + (x)] 


( 2 ) 


where the k's are constants and the w^ — are proportional to real 
solutions u y ±(x) (valid in the classical region) and v^—(x) (valid 
in the non-classical region). The constants of proportionality are 
determined so that the w^— remain well-behaved in passing through the 
turning point, x = 0 . 

For convenience we take the non-classical region to be to the left 

of the origin. The solution of Eq. (1) with c real was found by 

3 

Lommel to be: 

w v ±(x > 0) = A + u v ± = A + (£/Q) 1/2 J + (O 

(3) 

x -1 

where £ « / Q(x) dx, v » n + 2 , A + are constants, and the 
are Bessel functions of real argument. In the non-classical region 
this solution takes the real form 

w v ±(x < 0) = B + V v ±(x) = B + (U|/|Q|) 1/2 I +v (U|) (4) 

o 

where | £| = |Q(x)|dx , B + are constants, and I (z) = exp(-TTVi/2) J + ^(iz) 
are Bessel functions of imaginary argument. 
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Watson^ gives the following limiting forms of these Bessel 
functions : 

x -*■ 0 -■ 

( i 0±v 

= r(i + v) + 


<4 l5|) 


+v 


(5) 


r(i + 


V) 


+ ... 


X -*■ 00 


J + v ( ^ 


c^) 1 ^ 2 cos(£+ Vir/2 - tt/4) + 


( 6 ) 


x - 00 

I +V (U|) - (2tt | £| ) -1/2 { J^[l+OCj£j-)] + e I "l 5 l" ( I ± V)7Ti] + , 

(7) 

A derivation of the connection formulas, originally due to Jeffreys,"^ 
uses the limiting forms near the turning point. Expanding the w^~(x) 
of Eqs. (3) and (4) in powers of x , we obtain as leading terms 

w v + (x > 0) = A + (2v) 1/2 (cv) v x/r(l + v) + ... 
w v '(x > 0) 
w v + (x < 0) 

w "(x < 0) = B_(2v) 1/2 (cv)V(l - v) + ... 


= A_ (2v) 1/2 (cv)“ v /r(l - V) + ... 

= B + (2v) 1/2 (cv) v |x|/r(l+v) ■ + ... 


( 8 ) 
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A connection formula arises from the condition that these different 
solutions join smoothly at the turning point. This condition is satis- 
fied for w^ + (x > 0) and w^ + (x < 0) if B + = - A + and for (x > 0) 
and w^ (x < 0) if A = B__ . It is convenient to take the magnitudes 
of these constants as unity. Now the explicit expressions for the w^— 
in terms of the Bessel functions are 


w^ + (x > 0) = u v + (x) = C^/Q) 1 ^ 2 J +V (D 
(x > 0) = u v (x) = (C/Q) 1 ^ 2 J _ v (0 

w v + (x < o) = -v v + (x) = - (U|/|q|) 1/2 i +v (UI) 
w v (x < 0) = v v "(x) » (U|/|q|) 1/2 I_ v C | C l > 


O') 


(4') 


The asymptotic behavior of the 


and (7) : 



functions follows from Eqs. (6) 


w v + = O^rrQ) 1/ ^ 2 cos ( ^ - w/2 - tt/4) + ... 
= (-|itQ) 1 ^ 2 cos + w/2 - ir/4) + ... 


(9) 


x _ 00 


w v + - — ( 2 tt | Q | ) 1/2 { [1 + °(j||-)] + { e [ ^ ( 2 + V)7Ti] + ...} 

w v " - ( 2 tt | Q I ) -1/2 { J5 1 [x + 0 (y||-)] + { ~ V)lTi] + 


( 10 ) 
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We now wish to find linear combinations of these functions which 
are suitable for describing wave functions. If there is no turning 
point to the left of the origin, one combination must yield a decaying 
exponential in the non-classical region. The linear combination 
w v + + w v has the asymptotic behavior 

x - 00 

W V + + w v ” " 2 sin(VTT) (2 tt|q|)" 1/2 + ... (11) 

X <+ 00 

w v + + w^ =2 cos(w/2) (-^tQ) *^ 2 cos (5 - tt/4) + . (12) 

As usually stated, the connection formula in this case would have the 
form: 

sin(w/2) | Q | ly ^ 2 e ^ ^ Q ly/2 cos(£ - tt/4) o (13) 

The one-sided nature of this connection formula is meant to imply 

the following: if the expression on the left is a good asymptotic 

approximation to the true solution to the left of the turning point, 

then the expression on the right is a good asymptotic approximation to 

the right of the turning point. However, it is not possible to apply 

the converse, that is, to reverse the arrow. The reason for this bias 

in direction is not clear from Jeffreys 1 derivation (Jeffreys used a 

6 

double-head arrow), but is explained in Langer T s article where the 
case for V = 1/3 is treated. It is sufficient to state here that 
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when the arrow is reversed in Eq. (13) the increasing exponential term 
in w^— is introduced if the phase of the wave function is only 
slightly different from tt/ 4 . The connection formula problem is 
interesting from a mathematical viewpoint but is unimportant in practice.- 
In obvious analogy to the usual Airy functions we identify the 

(Xi^ as the linear combination w^ + + w^ and the 10i^ as the linear 

I “3/2 - 1/2 

combination w^ - w * Choosing = 2 csc(vtt) and = 2 

as the constants in Eq. (2), the following expressions are obtained: 
Oi v (-x) = \ csc(vrr) v 1/2 x 1/2 [I_ v (Q - 

&i v (x) = j CSC (vtt ) V 1/2 x 1/2 [J v (0 + J_ v (0] 

13i v (-x) = v 1/2 x 1/2 [I_ v (S) + I v (0] 

lgi v (x) = V 1/2 x 1/2 [J_ v (0 - J v (0] (14). 

1_ 

where, as defined before, £ = 2cV x . The multiplicative factors 

have been picked to make the Wronskian W of the generalized Airy 
functions the same as for the usual Airy functions 

W-Klli^(x) , I3i v (x)} = (Xi v (x) 10i v (x) - ^i^(x) ^ (X ° '' T_1 

(15) 

The generalized Airy functions, of course, obey the same type of rela- 
tions that the Ai and Bi functions do. In particular, the identity 


(*) - e ( " 1/2 -’™i) (k [xe 
V v 


( 16 ) 
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is operative- The series expansions for these generalized functions 
are readily obtained from the series representations of the Bessel 
functions. 

2 

Gordon has recently calculated the Ai and Bi functions to 
great accuracy using the method of generalized Gaussian quadrature. 

In order to apply this method to the calculation of the and 

^0i^ functions, it is necessary to express them in terms of the Bessel 
function of imaginary argument K^(z) 5=1 ttcscCvtt) [I ^(z) “ I^(z)]. 
Using the identity^ 


jr / \ 'I / \ 1/2 

K (u) = TT cos(VTT) u 


e ; o “172 — <17 > 

0 x /Z (jrfu) 


and Eq. (16), the pair of relations are obtained: 

n . t v _ 1 AwW -3/2 -1/2 1/2-1/ 4\) T p(u) , 

x) 2 cos(w)7T C X a /q 1+(u/^) ^U 

-a, / n , v -3/2 -i/2 “ p( u ) 

IOi v (-x) = cos (vrr)TT c x e I Q il, '' ( u /g ) du 


where 


P(u) = (2 /tt) 1/2 u 1/2 e" u K v (u) 


and c in the constant in Eq. (1), 


Unfortunately, Gordon's Eq. (A9) has a typographical error. It 


should read 


, . 0 -l/2 -3/2 -1/2 — x . . 

p(x) = 2 TT x e K 1 ^ 3 (x) 

Using this equation, his Eq. (A10) is obtained. Also his Eqs. (A14) 
and (A15) are wrong. However, Eqs. (A16) and (A17) are, correct. 
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In terms of , the generalized Airy functions for positive x 

are: 


0 -i v (x) = i /2 tT 1 csc(w/ 2 ) v 1/2 x 1,Z [K v (iO - K v (-iO] 
tli v (x) = 2 tt ~ 1 sin(w/ 2 ) v 1/2 x 1/2 [K v (iO + K v (-iO] 


( 20 ) 


The corresponding integral expressions are given by 


di v (x) = \ c 1/2 rr" 3/2 csc(vtt/2)cos (vtt) x 1/2 “ 1/4v /“ 




cos(g— tt/ 4) + (u/£)sin(£-Tr/ 4) ^ du 

1 + (u/5) 2 


'Sgi (x) = 2 c 1/2 tt 3/2 sin(vrr/ 2 )cos(w) x 1/2 l/4v / 
v 0 


( 21 ) 

where p(u) is the same as Eq. (19). 

The moments, y , of the positive function p(u) can be found 

g 

using the identity 





00 



X 


j-1 


-X 

e 


K^(x)dx =* tt 


- 1/2 


9 -j r(j+v)r(i-v) 
r(j+i/2) 


( 22 ) 
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Thus , 


y 


k 


/ x k p(x)dx 


r(k + 1/2 + v) rck + 1/2 - v) 
2 k k! 


( 23 ) 


Using the y , the integral expressions for the generalized Airy 

tC , 

functions can be approximated as sums of n points and n weights 

determined by a generalized Gaussian quadrature algorithm given else- 
9 

where. 
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